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AN EQUATION FOR VORTEX MOTION INCLUDING 
EFFECTS OF BUOYANCY AND SOURCES 
WITH APPLICATIONS TO TORNADOES 

By Robert C. Costen 
Langley Research Center 

SUMMARY 

A new equation is derived for the motion of vorticity in a general fluid, including 
the effects of viscosity, compressibility, nonhomogeneity, and nonconservative forces. 
This equation results from a kinematical revision in the derivation of Kelvin’s circulation 
theorem. It applies, in particular, to vortices which may not move with the fluid flow. 
The equation also admits fluid sources. The source term accounts for the influence of 
winds in one atmospheric layer on the motion of vortex tubes in an adjoining layer when a 
flux of fluid exists between the layers. 

A linearized form of the equation is used to obtain possible explanations for (1) why 
tornado cyclones move to the right of the mean tropospheric winds, and (2) why the 
observed revolution rate of the funnels of a twin tornado about their common center was 
retarded with respect to the fluid flow. This twin tornado occurred on April 11, 1965, 
near Elkhart, Indiana. 

It is shown that the rightward motion of tornado cyclones may be attributed on the 
basis of order -of -magnitude calculations to the upward flux of air from the core of the 
cyclone into the jet stream. It is also shown that a slowdown in the revolution rate of a 
twin tornado would result from (1) an attractive buoyancy force acting on each partially 
rarefied core in the presence of the radial pressure gradient of the other vortex or (2) an 
attractive magnetostatic force between axial dc electric currents of like sense in both 
cores. The retardation in revolution rate of the twin tornado due to buoyant attraction 
was estimated to be 5 percent. The electric current required for comparable retardation 
is on the order of 10® amperes, which appears to be much too large for atmospheric 
phenomena. According to the linearized theory, the maximum possible retardation for 
two buoyant vortices of like circulations and cross sections is 22 percent. This figure 
applies to hollow vortices revolving about each other in near contact. 



INTRODUCTION 


An accumulation of evidence exists that atmospheric vortices associated with severe 
storms do not necessarily move with the local fluid flow. According to reference 1, 
pp. 71-72, and reference 2, pp. 319-321, a storm cell which is destined to become severe 
moves with the mean tropospheric winds during the initial stages of formation. As the 
storm matures it acquires cyclonic rotation (being thereafter termed a tornado cylcone, 
see ref. 3, pp. 3-4) and commences to move at an angle of about 25° to the right of the 
mean tropospheric winds in the northern hemisphere. This is an example of atmospheric 
vorticity which does not move with the fluid. 

A second case of such deviate vortex motion is pointed up in this report. It con- 
cerns a twin tornado which occurred on April 11, 1965, near Elkhart, Indiana, and whose 
motion is documented in reference 4. Analysis of the data presented in this reference 
shows that the tornado pair apparently revolved about each other at a much slower rate 
than the fluid flow in which the tornadoes were immersed. 


The explanation for these motions would seem to lie in a theory embracing the 
motion of vortex tubes, such as the circulation theorem of Kelvin (ref. 5, pp. 35-37). 

Kelvin obtained a formula for the time rate of change of vorticity flux ^ ^ w • n dS 

through an arbitrary surface S of finite area which moves with the fluid, where 
C = curl v is vorticity and n is the unit normal on S. (See ref. 5, pp. 35-37, and 

ref. 6, pp. 150-151, 162. By Stokes’ theorem ^ J'J* Z? • n dS = v • d? which is the 

time rate of change of circulation about the closed circuit line l which bounds surface S 
and moves with the fluid.) When the fluid is inviscid and subjected to conservative forces 
only, and when the density is either a constant or a single -valued function of the pressure 

(autobarotropic fluid), this formula reduces to ^ 3> • n dS = 0, from which the con- 
clusion is reached that under these conditions, the vorticity moves with the fluid. The 
expression ^ ^ ZD • n dS was also considered by Helmholtz, V. Bjerknes, and others 

(as mentioned in refs. 7 and 8). V. Bjerknes evaluated the expression for the case of an 
inviscid nonbarotropic fluid and considered the effect of the Coriolis force (as discussed 

in ref. 9, pp. 233-261). The expression for ^ w • n dS in a viscous, compressible 

fluid is given in reference 10, pp. 51-52. 


But in all these treatments the surface of integration S (or bounding circuit 
loop l) is taken to move with the fluid. What is desired is to have S and l move with 
the vorticity, independent of the fluid, when these two motions are not the same. This 
distinction, being purely kinematical, is readily incorporated. The purpose of this report 
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is to modify the previous formulas for ^ w • n dS for application in the general 

case where surfaces of integration S may move in a manner different from that of the 
fluid, and to derive therefrom an unrestricted formula for the motion of the vortex tubes 
in a general fluid. A linearized form of this formula is then applied in rudimentary 
attempts to account for the two cases of atmospheric vortex motions cited. The formula 
may also apply to the motions of hurricanes and frontal cyclones, although these cases 
are not treated herein. A different approach to the movement of rotating storms is pre- 
sented in reference 11. 

All the usual terms of meteorological importance, such as the viscous terms and 
the Coriolis and centrifugal forces, are included in deriving the equation for vortex 
motion. In addition to these terms, fluid sources are included as a simple means of con- 
sidering the influence of winds in one atmospheric layer on the motion of vortices in an 
adjoining layer when a flux of fluid exists between the layers. 

The derivation of the magnetic and velocity fields at the ground plane near a tilted 
hydromagnetic vortex is given in appendix A. Appendix B is a description and tabulation 
of empirical data on the twin tornado of April 11, 1965, near Elkhart, Indiana. In appen- 
dix C the circulation of the twin tornado is deduced from its revolution rate. 

SYMBOLS AND NOTATION 
A cross-sectional area of vortex tube, m2 

a radius of funnel A of twin tornado J -2 at ground, m 

b ratio of circulations, Tg/r A 

c dimensionless parameter (eq. (42)) 

sin 2 x cos 2 ^^ + tan 2 x sin 2 i^ A (see eq. (70c)) 
electric field, V/m 
Lorentz force per unit length, N/m 
Lorentz force per unit volume, N/m 2 
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magnetic field, A/m 

axial electric current (figs. 11 and 12), A 

effective axial electric current (eqs. (71b) and (7 Id)), A 

electric current density, A/m2 

length, m 

circuit line of integration, m 

excess mass present in a vortex tube per unit length, kg/m 
unit normal 

center point of revolution of twin tornado J-2 at the ground (figs. 11 and 19) 
pressure tensor including viscous terms, N/m2 
scalar pressure, N/m^ 

pressure field which is single -valued function of density (eq. (23)), N/m^ 
deviation of pressure from single-valued function of density, (eq. (23)), N/m^ 
total fluid source strength, kg/sec 

mass of source fluid being added to a vortex tube per unit length per unit 
time (eq. (36)), kg/m -sec 

fictitious source (sink) strength per unit length in a vortex tube having con- 
vergent flow and updraft, kg/m-sec (see section entitled "Tornado model") 

electric charge, C 


source density; mass per unit volume and per unit time being added to a flow 
from a fluid source, kg/m^-sec 


radius from tornado touchdown point in ground plane, m 


center -to-center distance between touchdown spots of funnels A and B of 
twin tornado J -2, m 

distance of vortex A from center point of revolution P of twin vortices 
in ground plane (figs. 11 and 19), m 

distance of vortex B from center point of revolution P of twin vortices 
in ground plane (figs. 11 and 19), m 

surface area of integration, 

time, sec 

velocity of displacement of vorticity, m/sec 

arbitrary abstract velocity field utilized to describe the motion of regions of 
integration S and V, m/sec 

volume of integration, m3 

fluid velocity field, m/sec 

Cartesian coordinates, m 

angle shown in figure 16, rad 

dimensionless parameter (eq. (62c)) 

circulation of a fluid vortex, m^/ sec 

effective circulation of a fluid vortex (eqs. (70b) and (70e)), m^/sec 
small dimensionless parameter 
permittivity of vacuum, F/m 


angle depicted in figures 13, 14, and 15, rad 



K 


ratio of axial electric current to circulation for a hydromagnetic vortex 
(eq. (83)), A-sec/m 3 

A dimensionless parameter (eq. (61)) 

X,jn , v direction cosines 

I Lt 0 permeability of vacuum, H/m 

|, 77, £ Cartesian coordinates used in figures 13 and 14 

p mass density, kg/m 3 

p'(p) density field which is a single -valued function of pressure (eq. (19)), kg/m 3 

p" deviation of density field from single -valued function of pressure (eq. (19)), 

kg/m 3 

p e electric -charge density, C/m 3 

r vector shown in figure 16 

$ gravitational potential (including centrifugal force in the rotating frame of 

the earth), J/kg 

4> arbitrary scalar function, m 3 /sec 3 

<p azimuth angle in symmetrical coordinate system of figure 13, rad 

X tilt angle of vortex with respect to the vertical (figs. 4 and 8), rad 

if/ azimuth angle in ground plane of figures 15 to 18, rad 

^A’^B azimuth angles shown in figures 11 and 19, rad 

ij/k (I) retarded revolution rate of twin vortices with axial electric current I, 

rad/sec 

SI angular velocity of rotation of the earth, rad/sec 
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* w fluid v or ticity, curl v, sec _ l 

Superscripts: 

A funnel A of twin tornado J-2 (fig. 4) 

B funnel B of twin tornado J-2 (fig. 4) 

bl boundary layer at ground 

js jet stream 

0 unperturbed vector fields in perturbation expansion 

s source fluid 

1 first-order vector fields in perturbation expansion 

Subscripts: 

A funnel A of twin tornado J-2 (fig. 4) 

B funnel B of twin tornado J-2 (fig. 4) 

m moving point 

0 unperturbed scalar fields in perturbation expansion 

R , <p components in spherical coordinate system R ,cp,9 (figs. 13 and 14) 

R, components in cylindrical coordinate system (figs. 15 to 18) 

x,y,z Cartesian components of a vector 

C component in Cartesian coordinate system £,?], £ (figs. 13 and 14) 

1 first-order scalar fields in perturbation expansion 

± field component perpendicular to the ground or to the lower jet-stream 

boundary 



Mathematical notation: 


— vector (arrow is omitted for magnitude of a vector) 

( ) mean value over a time interval or space interval 

unit vector 

Dots over a symbol denote time derivatives. 

Designations: 

A funnel A of twin tornado J-2; touchdown point of funnel A 

B funnel B of twin tornado J -2 ; touchdown point of funnel B 

J-2 twin tornado of April 11, 1965, near Elkhart, Indiana 

DERIVATION OF GENERAL EQUATION FOR 
MOTION OF VORTICITY 

Fluid Sources 

Conservation-of-m ass equation jwith fluid sources in integral form.- The condensa- 
tion of water vapor in the atmosphere into rain or hail, which subsequently falls to the 
earth, is an example of naturally occurring sink flow (or negative source flow). The 
conservation-of-mass equation for the main fluid excluding the rain or hail is, in integral 
form (ref. 6, pp. 132-133), 

liT If dV + ^pv ■ fi CS = q dV (1) 

where the integrals extend over an arbitrary spatial volume V and its bounding sur- 
face S, with n taken as the outward unit normal on S. In this equation p is the den- 
sity of air and water vapor comprising the flow under consideration, v is its velocity 
field, and q is the rate at which mass density is being added to this flow from the source. 
(The source term q is negative for condensation, and positive for evaporation or 
sublimation.) 

Con servation-of -m omentum equation with fluid sources in integral form.- If the rain 
or hail, which will be termed the source fluid, may be described by a velocity field v s , 
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which is distinct from v, then momentum is also transferred from the source fluid to the 
main flow at the rate of qv s per unit volume. The conservation -of -momentum equation 
for the main flow becomes, in integral form (ref. 6, pp. 132-138, and ref. 9, pp. 233-236), 

iTJ ^(pv)dV + §[pv(v • n) -If • fTjdS = |qv s + grad $ + 2v X + ifjdV (2) 

where 

"p" pressure tensor (including viscous terms), N/m^ 


$ 


2p(v x n) 
f I 
T 


gravitational potential (including centrifugal force in the rotating frame of 
the earth), J/kg 

Coriolis force density, N/m^ 

angular velocity of the earth, rad/sec 

Lorentz force density, p e E + p 0 J x H, N/m^ 


In equation (2) the first integral represents the time rate of increase of momentum 
contained in spatial volume V, the second integral is the rate of efflux of momentum 
through the closed bounding surface S by convection and by the action of the surface 
tractions, and the third integral gives the rate at which momentum is being added to 
volume V by transfer from the source fluid, and by the action of the gravitational, 
Coriolis, and Lorentz forces. 

Conservation equations in differential form. - In differential form equation (1) is 


|£ + div pv = q (3) 

and equation (2) is 

2 

p ~ + v|£ + v div pv + p grad - p(v x u>) - div'p' = qv s + p grad $ + 2 pv x + F (4) 

where v^ = v • v and w = curl v (vorticity). When combined and divided by p, as is 
customarily done, these two equations give 

- v x co + grad - div^f + S.(v - v s ) - grad 4> - 2v x £2 - — = 0 (5a) 

0t & 2 p P v ' P 
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which for Cartesian coordinates may also be written 


^ + v • grad^ v = div p - ^(v - v s ) + grad $ + 2v x fi + ^ (5b) 

Vorticity Equation 

The vorticity equation is obtained by taking the curl of equation (5a) 

^ + curl j^(v x - ~ div "p" + j|(v - V s ) - 2v x fi - = 0 (6) 

An integral form of this equation is given by 

iX lir ' + XX ^ curl|j(v X u) - i divT + p - v s ) - 2v x fi - ^JdS = 0 (7) 

where the integrals are taken over an arbitrary spatial surface S of finite area, and n 
is the unit normal on S, as shown in figure 1. 



Figure 1. - Arbitrary surface S with unit normal n. 

Kinematic s 

Since the integrals in equation (7) are taken with time t fixed, the surface of inte- 
gration S may be regarded as moving and deforming in an arbitrary continuous manner. 
Not being a material surface, S need not move with the velocity v of the fluid medium. 
Instead, imagine surface S to be moving and deforming in accordance with some arbi- 
trary abstract velocity field U'(x,y,z,t), which will in general be different from the veloc- 
ity field v of the fluid. For such a moving surface of integration, the following kine- 
matical theorem of Helmholtz, as given in reference 7, pp. 130-132, applies: 
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( 8 ) 


A A • n dS = jj^ + U' div X - curl(tJ' x A)j • n dS 
where A(x,y,z,t) is an arbitrary vector field. In this equation the flux of A through 


A • n dS, is purely a function of time; 


the moving and deforming surface S, given by ^ 

hence, it is appropriate to take the total time derivative of this quantity as indicated in 
the first term of equation (8). 


General Kinematic Form of the Vorticity-Flux Equation 

If the vector field A in Helmholtz' theorem is taken to be the vorticity 3 = curl v, 
equation (8) becomes 


s 11 s • 8 dS = If [lr - curl ( 0, x 


n dS 


(9) 


Eliminating 83/at between equations (7) and (9) yields the general kinematic form of 
the vorticity -flux equation 


dt iX ^ “ dS + 


fT h • curl (u* - v) x 3 - i div p + ^(v - V s ) -2vXf2-^dS=0 


( 10 ) 


where the surface of integration S is convected with the abstract velocity field 
U'(x,y,z,t). The abstract velocity field U' and the fluid velocity field v which appear 
in this equation may be regarded as superposed, U* governing the motion of the surface 
of integration S, and v the motion of the fluid elements. These two fields of motion 
are (in concept) unrelated. Equation (10) reduces to the usual form of the vorticity -flux 
relation in the special case of setting U' = v. 


General Equation for Motion of Vorticity 

In relation (10) the abstract velocity U' with which the surface of integration S 
moves is continuous, but otherwise arbitrary. Now U' is constrained by the specifica- 
tion that S moves in such a manner that the flux of vorticity 3 through S is con- 
served; that is, 


Aj'jVSds-o 


(Ha) 


The field U’ when constrained in this manner is designated by the symbol U, where 
U satisfies the equation 

n • curl |(u - v) X 3 - i div"p + ^(v - v s ) -2vXfi-i|dS=0 (lib) 
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Thus, U may be regarded as describing the motion of the vortex lines which pass 
through S. 

Equation (12) applies to any surface of integration S having any chosen orientation, 
and the integrand is assumed to be continuous. This is possible only if 


curl 


(u - v) X co - — div p + 2(v - V s ) 


2v x O - -i-l = 0 


3 


or, equivalently, 




(u - v) x co - i div p + a(v - v s ) 

- 2 v x S2 - ^ = grad 0 


( 12 ) 


(13) 


where 0(x,y,z,t) is an arbitrary scalar field. 

Equation (13) is the desired general formula for the speed of displacement U of 
vorticity. Note that U is not, in general, determined uniquely because of the arbitrary 
function 0. Physically, this means that for general distributions of vorticity there are 
arbitrarily many flux-preserving motions. In the simple case when = Constant, any 
area-preserving motion of S perpendicular to Z o would also be flux conserving. But 
when a vortex tube is isolated, the displacement velocity U of the tube should be uniquely 
determined, except for displacements along the tube. (Also the term p grad <fi may 
absorb forces directed along oS which the term p(u - v) X cannot balance.) 

Substituting equation (13) back into the momentum equation (5a) gives 


a? 

at 


- U X co + grad(— — $ + 0] 


'(£-*♦*) 


= o 


or since in Cartesian coordinates 


2 

grad = ( v • grad) v + v xu 


then 


+ v • grad^ v = (u - v) x co + grad(4? 


(b) 


(14a) 


which shows that the Magnus force density p(u - v) x oo acts directly to accelerate the 
fluid elements. Comparison with equation (5b) shows that in equation (14a), all the non- 
conservative forces, including the viscous force, are lumped together into one term and 
represented by the Magnus force density. The vorticity equation (6) becomes upon sub- 
stitution of equation (13) 
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or 


^ - curl(u x 55) = 0 


^ + U • gradjco = (co • grad)u - u> div U (14b) 

which is analogous (with U replacing v) to relations first obtained by Helmholtz under 
restrictive conditions. (See ref. 7, pp. 129-130.) 

The pressure tensor Ip has been retained for generality to illustrate that a dis- 
placement velocity U may be assigned to vortex tubes even in the presence of viscous 
diffusion. For an inviscid fluid, equation (13) becomes 

(u - v) x w + i grad p + ^(v - v s ) - 2v x - ~ = grad (p (15) 

r r r 

where p is the scalar pressure. 

Alternative Forms of the Vortex-Motion Equation 
for an Inviscid Fluid 

The alternative formulas derived in this section are included here for completeness 
although they will not be used in subsequent sections. Because <p is arbitrary, any con- 
servative terms on the left-hand side of equation (15) (i.e., terms which have the form of 
a gradient of some scalar quantity) may be deleted from the left-hand side and absorbed 
in the arbitrary grad (p term. Therefore, since 

i grad p = grad^ - p grad^ (16) 

equation (15) for an inviscid fluid can be expressed as 

(tj - v) x 55 - p grad^j + j|(v - v s ) - Zv x £2 - i = grad c p (17) 

For an autobar otropic fluid, where pressure is a single-valued function of density (that is, 
P = p (p)), the term ~ grad p in equation (15) and p grad^j in equation (17) may be 

deleted by absorption into grad <p. If, in addition, all nonconservative terms should be 
omitted, these equations become 

(u - v) x 55 = grad (p (18) 


13 


from which it may be inferred by taking 0=0 that the components of U and v per- 
pendicular to co are equal, in agreement with Kelvin’s theorem that vorticity moves 
with the fluid under these conditions. 

The density may be written 

p = p’(p) + p” (19) 

where p' is a single-valued function of the pressure and p" is the deviation of p 
from p'(p). It follows that 


1 _ 1 p" 

P^pTpJ-j^ y 


( 20 ) 


and equations (15) and (17) may be written, respectively, 


(u - v) x oo - g ra -d P + jj(v - v s ) - 2' v x O - | = grad 0 


<t> 


where the terms 


(0 - v) x (5 + p grad p/ ^,( p )l + ^ (v - v s ) - 2v x 0 - | = grad 

[p’(p)] 


( 21 ) 


( 22 ) 


p’(p) 


grad p and p grad 


have been absorbed in the grad 0 


term. Alternatively, if the pressure is written in the form 


P = P'(p) + P" 


(23) 


where p’ is a single -valued function of the density and p" is the deviation of p 
from p’(p), equations (15) and (17) become, respectively, 


(u - v) x oo + i grad p" 

r 

+ *(v - v s ) -2 vXfl-^= grad 0 

(24) 

(u - v) x3 - p” grad^j 

^ + ■ l(v - v s ) - 2 v x - ^ = grad 0 

(25) 


From equation (13), (15), or (17) a speed of displacement U may be assigned to 
every element of vorticity in a fluid when the other relevant fields, such as pressure, 
density, and fluid velocity, are known. Evidently U frequently differs from v. These 
equations may prove to be useful in forecasting the movement of severe local storms and 
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tropical storms. This is because they show explicitly how various atmospheric fields 
affect the movement of vorticity, and vorticity is a characteristic property of such 
storms. Conversely, they may aid in determining the structure of such atmospheric 
vortices from observations of their motion. 


Linearization of Equation for Vortex Motion 

Perturbation series .- Analysis of equation (15) for the motion of vortex lines in an 
inviscid fluid is simplified by linearization as follows: 


v = v°(x,y,z,t) + ev 1 (x,y,z,t) 
p = P Q (x,y,z,t) + ep 1 (x,y,z,t) 
p = p Q (x,y,z,t) + ep 1 (x,y,z,t) 

4 > = 0 o (x,y,z,t) + e 0 1 (x,y,z,t) > 

w = ew 1 (x,y,z,t) 
q = eqjfoy^t) 

T = e? 1 (x,y,z,t) 


(26) 


where e is a small dimensionless parameter. Substitution of set (26) in equation (15) 
gives for the unperturbed potential flow 


j- Vp G - 2v° X n = V0 Q 

From equation (5a) (with viscous terms omitted), 


(27) 


Vd> = Vl<i> I _ 8v . 

9 o v V*o 2/ 8t 


(28) 


where €> Q is the potential for gravitational, centrifugal, and other conservative forces 
and Vq = v° • v°. The first-order perturbation equation is 

p Q (u - v°)x q^v 0 -V s ) + Vp x - 2p Q v 1 X O - 2pjV° X -’f 1 = P 0 ^4>i + pjVc^ (29) 
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Substituting equation (27) in equation (29) gives 

p G (u - v°) X co 1 + q^v 0 - v s ) + Vp x - ~ Vp Q - 2p o v! X O - T 1 = p (30) 

H o 

For the present purposes it is desirable to simplify equation (30) further by elimi- 
nating the perturbation pressure p^ and velocity v 1 fields. The term Vp \/p Q ma y 
be absorbed into the term provided the unperturbed density p Q is taken to be a 

constant. Elimination of v* requires that the Coriolis force be neglected. With these 
two restrictions equation (30) becomes 


p 0 (u - v°) x w 1 + q^v 0 - v s ) -p- Vp Q - I 1 = V(p 1 


( 31 ) 


In many cases the arbitrary term in equation (31) may be set equal to zero. But 

in general it may be needed in order to balance forces directed parallel to the vortex 
lines, which the Magnus term p Q (u - v°) x S* cannot balance. 


Physical interpre tation erf terms in linearized equation for vortex motion.- The 
terms of equation (31) all have the dimensions of force per unit volume. The first term 
P 0 (U - v°) x u>l represents the Magnus force density on vortex filaments which do not 
move with the unperturbed flow velocity v°. The second term q^(v° - v s ) represents 
the force density required to accelerate the source fluid from the source velocity v s 


to the flow velocity v°. The third term 


Pi 


Vp Q represents a "buoyancy" force which 


acts on density variations in the presence of the unperturbed pressure gradient. It corre- 
sponds, for example, to the buoyant lift on a balloon. The fourth term f * represents an 
extraneous, nonconservative force density, which in this study will be taken to be the 


magnetic force density T* = T* x B° where j 1 is a perturbation electric current den- 
sity, and B° is the applied magnetic field. 


T 1 


Equation (31) states that the vorticity moves in such a manner that the Magnus 
force, source -force, buoyant -force, and magnetic -force densities vanish (or are equal to 
the gradient of some scalar field). Only the Magnus term p Q (u - v°) x £>■*- contains the 
velocity of displacement U of the vorticity. Incidentally, in the derivation of equa- 
tion (31) no attempt has been made to follow the motion of any quantity except the vorticity. 
That is, source strength is not necessarily conserved and density irregularities are not 
necessarily followed by the motion tT. 


Linearized equatio n for the motion of isolated vortex tubes. - Up to this point the 
vorticity has been considered to be distributed generally throughout the flow. The case 
of an isolated vortex tube of circulation r will now be considered as a prelude to the 
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treatment of atmospheric vortices. For application of equation (31) to an isolated vortex 
tube, it is convenient to integrate over the volume V of the tube included between two 
surfaces S' and S" which are everywhere normal to the vorticity in the tube, as shown 
in figure 2. 



Figure 2.- Isolated vortex tube with cross-sectional 
surfaces S' and S" and included volume V 
shown. 

Upon integration, equation (31) becomes 


OTk(H - v°) X w 1 + q,(vO - v s ) - ^ Vp Q - T 1 - V0 1 


dV = 0 


(32) 


This volume integral may be written as an integral on a cross-sectional surface S and 
along a line l parallel to the vorticity: 


J d l dS p c (u - v°) + qi (v° - v s ) - ^ Vp Q - T 1 - V0 


= 0 


(33) 


Since the integrand is assumed to be continuous throughout V, mean values may be 
assigned to the quantities (u -v°), ( v° - v s ), Vp Q , and on surface S so that 

equation (33) is termwise identical with the following equation: 

J di jp 0 <l! - v°> x JJ oJl dS + <v° - v s > jj q x dS - P x dS 


-H Tl dS " 



(34) 


where the brackets ( ) denote mean values on the cross-sectional surface S. This 
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equation may also be written 

§ p Q (u - v°> x f 1 + (v° - v s )Q 2 - 

where 


M t - F 1 


< v *i) 


dz = 0 


circulation of the vortex tube 

mass of source fluid being added to the vortex tube per unit length 
excess mass present in the tube per unit length 
F* Lorentz force on the tube per unit length 

A cross-sectional area of the tube (taken perpendicular to 


Ql 

Mi 


(35) 


The remaining integration in equation (35) is over an arbitrary length l of the 
vortex tube, and the integrand is assumed to be a continuous function of l ; hence, the 
integrand itself must vanish for all values of l : 


P D <U - v°> X f 1 + <v° - v s >Q 1 - ^ p P °^ M x - F 1 = ( V<^ 


(36) 


The arbitrary term A has been retained for the purpose of balancing possible 

axial forces which cannot be balanced by the Magnus force p Q ( U - v°) x ?1. This first- 
order perturbation equation is applicable to an isolated weak vortex tube in an inviscid 
fluid without Coriolis forces. The unperturbed flow must be irrotational and free of 
sources, have constant density p Q , and must satisfy the equation 



(37) 


where <f> 0 is the potential for gravitational, centrifugal, and other conservative forces. 

In subsequent sections where equation (36) is applied, the subscript 1 (or super- 
script 1) will be dropped from the perturbation fields T*, Qp Mp and F*. The sub- 
scipt o (or superscript o) will be retained on the unperturbed fields, however. Any 
confusion which might arise from this practice can be resolved by referring back to 
equation (36). 
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APPLICATIONS OF LINEARIZED EQUATION FOR VORTEX MOTION 
TO TORNADOES AND TORNADO CYCLONES 


Rudimentary Approach To Be Used in the Following Applications 

In the following sections equation (36) will be applied to tornadoes and tornado 
cyclones as if they were weak perturbations on an unperturbed potential flow. All non- 
linear effects will clearly be missing in this rudimentary treatment. Moreover, only 
gross features of these storms, such as circulation, updraft, and core density, will be 
considered together with some effects of the boundary layer at the ground and of the jet 
stream near the tropopause. 


This simplified version of a severe local storm which is producing a tornado may 
be described briefly as follows: A tornado cyclone, or parent storm, is regarded as a 

vertical vortex which significantly affects the flow in a region about 16 km (10 miles) in 
diameter. The tornado cyclone has a core that is several miles in diameter and a 
strong updraft within the core. The cyclone terminates near the tropopause and the 
updraft from its core is ejected into the jet stream and carried eastward. A tornado 
produced by this storm is pendent from the cloud base at an altitude on the order of 
0.3 km (1000 ft) and may extend to the ground. The circulation of the tornado is about 
1 order of magnitude less than that of the tornado cyclone, and its core diameter is from 
1 to 2 orders of magnitude less. An updraft exists in the tornado core, and retarded 
fluid from the viscous boundary layer at the ground is drawn up into the tornado. This 
interaction with the boundary layer at the ground is considered capable of affecting the 
speed of displacement of the lower portion of the tornado - but not of the tornado 
cyclone. The motion of the tornado cyclone, however, may be influenced by its inter- 
action with the jet stream above. Since the boundary layer at the ground acts as a source, 
and the jet stream acts as a sink, both these effects depend upon the source term 
<V° - V s ) Q x in equation (36), and these are the first cases which will be treated. 


Because of the high speed of the circulating wind in a tornado, the tornado core may 

Oo> 


be rarefied so that its motion is influenced by the buoyancy term 


Mj in equa- 


tion (36). Also, the updrafts in a tornado cyclone indicate that the core of the cyclone 
may be warmer than the surrounding air and therefore buoyant. If the vortices should be 
tilted from the vertical, one effect of their buoyancy would be the presence of upward lift, 
which could affect their motion. This effect, although potentially important, will be 
omitted in the case of the twin tornado because the algebra becomes cumbersome. It 
will also be omitted from the treatment of tornado cyclones because it appears that in this 


1 


For a more complete discussion of severe local storms, see reference 2. 
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case the term Vp^ in equation (30) cannot be absorbed into the p 0 V</>^ term, as done 
throughout this report, nor can it be ignored. The case of buoyancy which will be treated 
with equation (36) is that of buoyant attraction between the rarefied cores of a twin 
tornado. 

The applications of equation (36) are concluded with a hypothetical example which 
depends upon the Lorentz term. The case treated is that of two parallel vortices which 
are subject to magnetic attraction by virtue of dc axial electric currents of like sense in 
the vortex cores. 

Effect of the Boundary Layer at the Ground on Tornado Motion 

Description and analysis .- For this analysis the following simplified form of equa- 
tion (36) is used: 

p 0 <U - v°> x F + <v° - v s >Q = 0 (38) 

where only the Magnus term and source term have been retained, and where the sub- 
script 1 (or superscript 1) has been dropped from the perturbation fields T and Q. 

Near the ground the flow undergoes a rapid transition from the free -stream flow to the 
retarded flow in the viscous boundary layer. Ordinarily, the boundary layer does not 
significantly affect the free -stream flow. But in the event of a tornado, the retarded 
fluid in the boundary layer is sucked up into the free stream as shown in figure 3(a). 

This constitutes a fluid source at the boundary of the free stream, which is coincident 
with the tornado and can affect its motion. 

The upward flux of fluid from the boundary layer into the tornado core occurs 
essentially as follows: The pressure p in the boundary layer is the same as in the 
free stream, although the flow is considerably retarded. Therefore the radial pressure 
gradient within the free-stream vortex also exists in the boundary layer below. The 
boundary -layer fluid moves radially inward under the influence of this pressure gradient 
and ultimately ascends into the low-pressure tornado core. Although the magnitude of 
this flux from the boundary layer to the free stream may not be accurately known, there 
is ample evidence of its existence (ref. 12, p. 675, and ref. 3, p. 7). 

This upward flux of fluid from the boundary layer appears as a source to the free 
stream. The total magnitude of this source, in kg/sec, is 

ZL = P 0 Av i (39) 

where A is the area of the tornado core at the free-stream boundary, and v_j_ is the 
upward component of velocity at this boundary. The fluid source is located (concentrated) 
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on the lower boundary of the free stream. From equation (38), therefore, a strong dis- 
placement, or kink, would be expected to develop in the tornado at this boundary. But 
observation of tornadoes indicates that this does not occur and that the tornado maintains 


v 1 - 



-Vortex tube 


_ v bl iV 

777777/ ////// 777 7777 

(a) Vortex tube drawing fluid up from the boundary layer. 



00 Top view of displacement velocity U of vortex tube relative to the 
free-stream velocity "v° and the boundary -layer velocity ~v^l. 

Figure Effect on vortex motion of the boundary layer at the ground. 


a smooth contour, as if the source were distributed over a finite length of the vortex tube 
instead of being concentrated at its lower extremity. It may be that there is a stabilizing 
effect in the tornado which acts to prevent strong curvature of the core, and in so doing, 
effectively distributes the concentrated source. 

Such a stabilizing effect is assumed to exist and the source strength is consid- 
ered to be distributed over a vertical length L of the lower part of the tornado, whence 


Q - 


Pq A v 1 

L 


(40) 
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and equation (38) becomes 


(tJ - v°) x F + (v° - v bl )^p = 0 (41) 

where v b * is the unperturbed velocity in the boundary layer (source fluid), and where 
the brackets have been dropped by taking mean values over the vortex cross section to be 
approximated by values at the core center. If 

vbl = cv° (42) 


with 0 < c < 1, then 


Av , 


(u - v°) x f + v°(l - c) — i = 0 


(43) 


In Cartesian component form, with circulation r taken along the z-axis and v° taken 
in the x-direction, equation (43) becomes 


Av, 

U y T+ v°(l - c)-^= 0 


(44a) 


(u x -v°)r=o 


(44b) 


whence 


U x = V o 


Av, 

u y = " v °(i ~ c ) ~rt 


(45a) 

(45b) 


It is clear from equations (45) that the effect of the boundary layer at the ground is 
to make the base of the funnel move to the right of the local wind v°, as shown in fig- 
ure 3(b). In the case of the twin tornado referred to in the Introduction and in appendix B, 
this implies that the lower portions of the two funnels should tend to separate as they 
revolve about each other. According to figure 4, which was prepared by Fujita (ref. 4), 
they did apparently separate from a center-to-center distance of 118.9 m (390 ft) to 
147.1 m (482.6 ft) in a period of 52 sec during which time the pair revolved about their 
common center through an angle of about 140°. 
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(a) Damage path and dehris marks near the Midway Trailer Court indicated 
by letter M. Six pictures were taken in the direction of the arrows. 



(*) Sequence depicting the initial occurrence of funnel A, the subsequent 
occurrence and growth of funnel B at the expense of funnel A, and the 
revolution of funnels A and B about their common center. 



Figure 4 .- Twin tornado J-2 path and configuration, as drawn by Fujita (ref. 4 ) 
from an aerial survey and from the six photographs taken by Paul Huffman. 
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Order -of -magnitud e check on boundary-layer effect. - In order to check whether the 
observed rate of separation of the twin tornado could possibly have been due to boundary - 
layer interaction, an order -of -magnitude calculation for this effect may be made as fol- 
lows: For the tornado pair the ratio of radial speed of displacement -Uy to tangential 
speed of displacement U x is on the order of -Uy/\J X ~ 0.1. From equations (45), 

-U v Av, 

_JL = (1 _ 0.1 

U x UL 


or 


v 


1 " 


TTL 

A(1 - c) 


X 


0.1 


(46a) 


A value for A, the area of the core at the base of one of the funnels, may be obtained 
from figure 5 

A = 6 x 10 3 m 2 


where the scale length used in this estimate is the center -to -center distance as given in 
table I, location 4. The circulation T of each funnel, when approximately equal in 
strength at location 4, is estimated in appendix C as 

r * 5 x 10 3 m 2 /sec 

Taking the mean velocity in the boundary layer equal to one -half the free -stream velocity 
gives, from equation (42), 

c = 0.5 


When these values for A, T, and c are substituted, equation (46a) becomes 

▼j. = ^ (46b) 

which is a relation between the updraft velocity v^ in the core at the boundary layer and 
the length L at the base of each funnel, which is assumed to be affected by the boundary - 
layer interaction. Taking L = 60 m (196 ft) gives v_^ = 10 m/sec (22.4 mph) for the 
updraft velocity at the boundary layer, a value which appears to be of reasonable magni- 
tude. It is concluded that the increasing separation observed for the twin tornado may 
have been due to its interaction with the boundary layer at the ground. 
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Figure 5*“ Photograph of twin tornado J-2 taken at 5:30 p.m. CST on April 11, 

1965 ; number 4 of a series of six pictures taken by Paul Huffman, staff photog- 
rapher of the Elkhart Truth. Left funnel is shown over the Midway Trailer 
Court (fig. 20). Reproduced by permission. 


Effect of the Jet Stream on Motion of Tornado Cyclones 

T he jet stream as a sink fluid . - In this section, equation (36) is again considered in 
the simplified form of equation (38) 

p Q <U - v°> x f + <v° - V s ) Q = 0 

where all terms except the Magnus term and source term are neglected. When a thunder- 
storm matures into a severe storm, which may generate tornadoes and hail, the storm 
cell begins to rotate cyclonically (in the northern hemisphere) according to observations 
reported in reference 3, pp. 3-4. This vortex is termed a "tornado cyclone." In addition 
to circulation, there is a strong axial flow of rising air within the vortex which frequently 
penetrates the tropopause. In the vicinity of the tropopause at midlatitudes is a narrow 
band of high-speed flow called the jet stream (ref. 13, p. 553). Much of the fluid in the 
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updraft of the storm is swept eastward by the jet stream as shown in figure 6, and forms 
the broad anvil canopy which is characteristic of these storms (ref. 1, pp. 67-68, and 
ref. 2, pp. 185-195). For the tropospheric flow below, therefore, the jet stream consti- 
tutes a sink fluid. 


V J' S 



Vortex tube- 1 

Figure 6.- Vortex tube ejecting fluid into the jet stream. 

The efflux of fluid from the tornado cyclone into the jet stream will affect the 
motion of the tornado cyclone in the troposphere below the jet stream. In particular, an 
attempt will be made to show on the basis of equation (38) that the storm will tend to move 
to the right of the mean winds in the troposphere in agreement with observations reported 
in reference 1, pp. 71-72, and reference 2, pp. 319-321. 

The winds in the troposphere below the jet stream vary considerably in direction 
with altitude, as described in reference 2, pp. 319-336. Moreover, the actual location of 
the fluid sink is near the upper extremity of the storm, where the cyclone core intersects 
the lower boundary of the jet stream. Despite the wind variation and sink concentration, 
the storms are observed to move, more or less, as a unit, although the core of the storm 
may be tilted from the vertical by as much as 20° (ref. 1, p. 69). Possibly the storm 
structure resists deformation of the core of the cyclone and in so doing effectively spreads 
out localized effects over its entire vertical extent. If this is the case, it would seem 
appropriate to distribute the sink strength uniformly over the vertical extent L of the 
storm and also to define a mean wind velocity v° for the troposphere below the jet 
stream. This is the view which is adopted in this section. The vortex tube is also 
assumed to flare out horizontally upon penetrating the jet stream so that the tube as such 
is confined to the flow region below the jet stream. 

The total strength of the sink is p Q Avj_, which in terms of total source strength 
in kg/sec is 
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(47) 


Z = -p 0 Av i 


where A is the cross-sectional area of the rising column of fluid at the lower jet- 
stream boundary and is its upward velocity at the boundary. When this source 
strength is distributed over the height L of the vortex, the source strength per unit 
length is 


Q = - 


Pq A v 1 

L 


(48) 


and equation (38) becomes 


(u - V°) X f - (v° 



(49) 


where v^ s is the fluid velocity in the jet stream (sink fluid) and v° is the mean fluid 
velocity in the troposphere below. The brackets are eliminated in equation (49) by taking 
the mean values over the vortex-core cross section to be approximated by values at the 
core center. If v° is taken along the x-axis of a Cartesian coordinate system and T 
along the z-axis, equation (49) becomes in component form 


u y r - (v° 




(50a) 





whence 


. Av . 

U x = v° + vi s 

x y n 


(50b) 


(51a) 


U 


y 





(51b) 


Sample solutions for the motion of tornado cyclones .- The dimensionless constant 
Av^^rL in equations (51) may be estimated by using the following typical values: 

Circulation of the tornado cyclone (ref. 1, p. 76), 

T = 10 5 m^/sec (52a) 
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Vertical extent of storm (ref. 2, pp. 176-177), 


L = 12.2 km (40 000 ft) 

Upward velocity of fluid in the cyclone core (ref. 12, p. 675), 
v ± = 152 m/sec (500ft/sec) 

Area of core, assumed circular with radius \/2 km (ref. 1, p. 74), 
A = 6.28 x 10 6 m 2 



(52b) 


(52c) 


(52d) 


(53) 


In order to determine the storm displacement velocity (T, two other quantities must 
be specified in equations (51). These are the ratio vj s /v° of the wind magnitude in the 
jet stream vl s to that below v°, and the included angle between v° and v3 s . (In the 

vector diagrams of figure 7, the vortex is viewed from above and v° is always taken to 

v js ^ 

lie in the same direction.) In figure 7(a), the wind-speed ratio is - — = £ and vector 

v° 2 

diagrams are sketched for various included angles between v° and VI s in 15° incre- 

v js 2 

ments. Similar diagrams are presented in figure 7(b) for — — = t-. Note that the speed- 

— v° 1 

of -displacement vector U is usually found to lie to the right of both v° and vl s in 
agreement with observation. The general condition that U lie to the right of v° is 
that (v3 s - V°) • v° > 0. 

Effect of Core Buoyancy on Revolution Rate of a Twin Tornado 

Occurrence of t win tornado J-2.- On April 11, 1965, a twin tornado (fig. 5) occurred 
near Elkhart, Indiana. Full details of this tornado (designated J-2) are presented in ref- 
erence 4, and an abbreviated account is given in appendix B of this paper. The funnels 
were parallel, about 122 m (400 ft) apart (center -to -center distance), and were tilted at 
an angle of about 29° from the vertical. Their circulations were of like sense (cyclonic) 
and they revolved about each other while slightly increasing the distance of separation. 

The relative motion of this tornado pair is an interesting case to which equation (13) 
should apply. Rather than attempt this at present, the linearized form, equation (36), will 
be applied, despite its limitation to weak vortices in an inviscid fluid without Coriolis 
forces. 
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Tornado model .- Figure 8 depicts a physical model for the twin tornado to which 
equation (36) may be applied. The vorticity is assumed to be largely contained in two 
tubes which are finite in cross section, semi-infinite in length, and terminate at the 
ground. The tubes are parallel, but tilted with respect to the vertical, and do not neces- 
sarily have the same circulation strength. Updrafts are present within the cores as an 
outlet for the convergent flow in the boundary layer at the ground, which was mentioned 
in the section entitled "Effect of the Boundary Layer at the Ground on Tornado Motion." 
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Figure 8.- Model used for twin tornado: tilted vortex tubes with 

updraft and influx. 

One model for tornadoes (refs. 14 and 15) assumes that such a convergent flow 
exists not only within the boundary layer at the ground, but at every point along the ver- 
tical extent of the funnel in the free stream. This convergent flow is conceived as main- 
taining the concentration of vorticity in the funnels. Thus besides circulation r, each 
core has convergence Q*, in kg/m-sec, which represents the rate at which convergent 
fluid is entering the core per unit length in the free stream. The asterisk notation is 
needed because within the core this convergent flow of fluid turns upward in a rising 
column which remains a part of a total free -stream flow field. If the converging fluid 
were actually withdrawn from the flow (as by entering a pipe), the asterisk notation would 
be unnecessary. 

This circulation-convergence -updraft model will be adapted herein for tornadoes. 
Thus the cores in figure 8 are assigned individual circulations and Tg and con- 
vergences and Qg. For a given location of the tubes the velocity field v outside 

the cores may be solved for, assuming that the external fluid is inviscid, incompressible, 
and irrotational, and that it does not penetrate the ground. The velocity field thus deter- 
mined should be a good approximation to the real flow except within the cores and within 
the boundary layer at the ground. The effect of the boundary layer at the ground, which 
is to make the cores move to the right of the local flow (and thus tend to separate) has 
already been discussed in the section entitled "Effect of the Boundary Layer at the Ground 
on Tornado Motion." 

Absence of source force density despite possible convergent flow. - The motion of 
the cores may be determined from equation (36). The source term (v° - v s )Q in this 
equation would apply if and Qg were true sinks. However, they are fictitious in 
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that the converging fluid is not withdrawn from the free-stream flow, but is simply 
deflected upward. This rising column does not constitute a sink except at a point where 
it penetrates through a surface of discontinuity into a different flow region. Any such 
penetration point is assumed to occur only at altitudes which are large compared with 
the separation distance of the two tornadoes, so that the relative motion of the two funnels 
near the ground is unaffected by it. Therefore, despite convergence in the tornado model, 

the source term (v° - v s ) Q in equation (36) is properly equated to zero for this study 
of the mutual interaction between the funnels. 


Buoyancy terms .- Now if the Lorentz force density F in equation (36) is also not 

< V P 0 > 


considered, the buoyancy term 
pair: 


M remains to influence the motion of the vortex 


p Q <U - v°> x f - M = A < Vcf > > 


(54) 


The possibility exists that the cores of the tornado pair are partially rarefied. This will 
occur if the circulating winds approach the speed of sound near the cores, as illustrated 
in the compressible -vortex example presented in reference 10, pp. 158-159. Thus the 
excess mass M per unit length is negative, and each vortex is buoyant in the radial 
pressure gradient of the flow field of the other. Consequently, the buoyancy force acting 
upon the pair should be attractive and tend to slow down their revolution rate. 

Ef fect of tilt neglected .- Although the convergences and Qg do not appear 

explicitly in equation (54), they still influence the velocity fields and pressure gradients 
in this equation and thus complicate the solution. The problem of two buoyant vortices 
with influx is simplified if the pair are taken to be alined with the vertical and to be iden- 
tical in influx and circulation, as illustrated in figure 9 where = Tg = T and 
Qa. = = Q* The equal influx and circulation strengths correspond to the state of the 

tornado pair at location 4 of figure 4. Neglecting the tilt angle appears to be warranted 
in the following order -of -magnitude calculation to determine whether buoyant attraction 
could cause signficant retardation in the revolution rate of twin vortices with influx. In 
appendix C it is estimated that the revolution rate of the twin tornado was retarded by at 
least 50 percent with respect to the fluid. 

Adaptation of linearized equation for vortex motion to t win -vortex model . - In order 
to calculate the motion U" of vortex A, the unperturbed flow, represented by v° 
and p Q in equation (54), is taken to be the flow field due to vortex B. Similarly, the 
motion of vortex B is determined by considering the unperturbed flow to be that due to 
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Figure 9*- Vertical vortex tubes having identical circulation 


and influx with centers at 


(x A ,y A ) and ( V y B ). 


vortex A. For either vortex equations (54) and (37) apply: 


P 0 <u 


V°) 


x r 



= A <V0> 


(55) 


where v^ = v° • v°. In the Cartesian coordinate system of figure 9 where the z-axis is 
vertical, the x- and y-components of equation (55) are, for vertical vortices, 


Po r z 

M 


/ \ 9v 

fa -v 0 \ + v 0 -^+v° 

V y y) * ax y 


irL 

8x 


at 


• = 0 


(56a) 


P o T z 

M 




8v° 8v° 

v° —X + — 21 = 

y 9y at 


(56b) 


where the brackets have been eliminated by taking the mean values over the vortex cores 
to be approximated by values at the core centers. The A (v<£>) term was retained in 
equation (55) only for the purpose of balancing possible axial forces. For vertical vor- 
tices this term need have only a z -component and therefore does not appear in equa- 
tions (56). 

Suppose now that the vortex giving rise to the unperturbed flow is located at the 
moving point x m (t),y m (t). The velocity components for the unperturbed flow are then 
given by 
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( 57 ) 


o _ K ( x ~ ^Sn) G (y-y m ) 
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where 


K 


2tt 


G = 


_r 

2v 


| 7 - 7 m|=/( x -x m ) 2 + (y-y m ) 2 


These velocity fields and their x-, y-, and t -derivatives, if evaluated on the line 
y - y m = 0, take the form 


v° = - 

K 

(58a) 

x ( 

x - x m) 

<_ 

K 

(58b) 

0X 

( x - x m) 2 

K 

G 

(58c) 

8 , 

( x - x m) 2 


*m K , y m G 

(58d) 

at 

( x " x m) 2 ( x " x m) 2 

v° - - 

G 

(58e) 

y ( 

: x ’ x n.) 

9v y _ 
ax 

G 

( x - x m) 2 

(58f) 

!ZL 

K 

(59g) 

»y 
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Ill 


mi mi 


!zL 

at 


( x - x m ) 2 ( x - Xm)' 


where dots denote time derivatives, and where K and G are taken to be constants. 

Now take the two moving vortices of figure 9 to be momentarily located on the 
x-axis at equal distances Rg/2 from the origin, as viewed from above in figure 10. 


a (core radius) 


(-M 


9 

/Rc 
\ 2 ’ 


Figure 10.- View of figure 9 from above when both vertical vortices are 
momentarily located on the x-axis at equal distances from the origin. 


Substituting equations (58) into equations (56) for the motion U A of vortex A in the fields 
due to vortex B gives 

K*_S* (59a) 

M Rc) ? R c 3 R c R 2 




x R, 


d 2 p2 


where the identifications x m = and y m = have been made. The symmetry of 
the problem dictates that U B = -U A . Hence, equations (59) become 


27rp o G / A _ G_\ K 2 G 2 U x K U y G _ 
M \ y Rcj R? Rc 3 R 2 Rp 
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(60b) 


27 T Po G / tt A K 


M 




U A G U, A K 


*4 


L =0 


Ri 


The excess mass M in the partially rarefied cores is negative and may be 
expressed as 


M = - A-rra^ 


P o 


(61) 


where a is the core radius and A, which has the range 0 ^ A ^ 1, is a dimensionless 
parameter which describes the degree of rarefaction. Substituting equation (61) into 
equations (60) and collecting terms gives the following equations for the displacement 
velocity U A of vortex A when it is located on the x-axis at a distance of R,, meters 
from the vortex B (also on the x-axis): 


uA (1 + £ ) + u AM-|L + gA + g(Uo 


R c GR C 


Rn 


(62a) 


U~(l - j3) + U. 


A K£ K_ _ 

y g " r„ 


0 


(62b) 


where 


0 = 


Aa^ 


2r: 


(62c) 


— A 

Solving these equations for U gives 





1 



(63a) 


(63b) 
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Result (63a) shows that one effect of buoyancy is to speed up the rate of approach of two 
vortices with convergent flow (Q* negative). It will be shown from result (63b) that 
another effect of buoyancy is to slow down the revolution rate of vortices of like sense. 

Correlation of t he twin-vortexjnodel with data from the twin tornado.- According 
to results (63) for the model of twin vortices with convergence, the two vortices spiral 
inward toward each other, and buoyancy serves to increase the rate of approach. But 
instead of spiraling inward, the two funnels of the twin tornado actually increased their 
distance of separation from 118.9 m (390 ft) to 147.1 m (482.6 ft) according to figure 4(b). 
Therefore, the convergence Q* of the twin tornado may be regarded as zero, or as dom- 
inated by the boundary -layer interaction which tends to make the vortices separate. This 
simplified comparison, which neglects tilt, casts some doubt on the existence of significant 
convergence in tornadoes except within the boundary layer at the ground. 

The revolution rates predicted by the present model and those observed for the 
twin tornado will now be compared. If Q*^/r^ is neglected, as seems permissible 
from the preceding analysis, result (63b) becomes 


U. 


r/2^R c 



(64) 


With both vortices momentarily located on the x-axis, equation (64) gives the ratio of the 
revolution rate of a buoyant vortex pair to what it would be if the vortices moved with the 
fluid. This expression will now be evaluated for the time when the tornado was in loca- 
tion 4 of figure 4, at which time both funnels appear to have had approximately equal 

strengths (in agreement with the assumptions of this section) according to the sketches in 

a 2 

figure 4(c). From the photograph of the pair at location 4 in figure 5, — « 0.1 at the 

ground. If A is taken to be 0.5, which corresponds roughly to sonic speeds at each core 
periphery, buoyant attraction is found to retard the revolution rate of twin vortices by 
about 5 percent. 

According to an estimate made in appendix C, the revolution rate of the twin tornado 
was retarded by at least 50 percent with respect to the fluid flow. The effect of buoyant 
attraction, although tending in the right direction, is not sufficiently strong to account for 
this much retardation in the linearized theory. 

The maximum possible value of a^R^ for two vortices of equal core size is 1/4, 
which corresponds to the two vortex cores in actual contact. The retardation in this case 
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for A = 0.5 is then about 12 percent. Totally evacuated cores (A = 1) or cores filled 
with hydrogen gas (A * 1) would raise the maximum retardation to about 22 percent. 

Influence of Axial Electric Current on Revolution 
Rate of Twin Vortices 

Speculative hydromagnet ic -vortex hypothesis . - In this section equation (36) will be 
taken in the simplified form 

p Q <U - v°> X f - F = 0 (65) 

The Lorentz force term F in this equation is capable of exerting an attractive force on 
the funnels of the twin tornado which could retard their revolution rate to any degree, 
stop it altogether, or even reverse it. The possible existence of such electromagnetic 
forces in tornadoes is highly speculative, but it does provide an instructive application of 
equation (36). For a bibliography of other electrical theories concerning tornadoes, see 
references 16 and 17. 

It is well known from electrodynamics that two parallel wires will exert an attrac- 
tive force on each other if (a) they are oppositely charged or (b) they are conducting elec- 
tric currents of like sense. Since the two funnels of the twin tornado were in contact 
with the ground, it seems unlikely that opposite static electric charges on the two funnels 
could be maintained. But axial electric currents of like sense could conceivably have 
existed in the two funnels, and this is the case that will be treated in this section. 

Fluid vortices with strong axial dc electric currents are termed hydromagnetic 
vortices. They have been discussed by Busemann (refs. 18 and 19) for the somewhat 
different case of a fluid with uniformly high electrical conductivity. In the analysis to 
follow high conductivity is assumed to persist only in the tornado cores. 

Hydr omagnetic -vortex model .- Each vortex tube of the twin-tornado model (fig. 8) 
is characterized by coincidence of the fluid and electromagnetic states depicted in fig- 
ure 11. This sketch is a view of the ground plane in the proximity of the touchdown 
points A and B of the tornado funnels. The broken lines represent projections in the 
ground plane of the tilted axes of tornadoes A and B. 

Equation for vortex motion with axial electric current i n a magnetic field .- The 
Lorentz force per unit length on an electric current filament in the presence of a mag- 
netic field is given by 

F = ju 0 T x H (66) 
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where I is the current in amperes, H is the magnetic field in amperes/meter, and 
is the magnetic permeability of vacuum in henries/meter. Substituting equation (66) 
into equation (65) gives 

p Q <U - V°> X f + fX 0 <H°> x T = 0 (67) 

where in this linearized equation, T is viewed as a perturbation current and H° as 
the unperturbed magnetic field. 



4 > 

B 

(a) Velocity fields in ground plane due to tilted vortex tubes 
of circulation pA 

and r B . 





B 


(b) Magnetic fields in ground plane due to tilted current 
filament s I A and I B . 


Figure 11.- Duplicate view from above of the ground plane , showing touchdown 
points of tilted hydromagneti c vortices A and B, with projections of the 
vortex axes indicated by broken lines. 
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As in the section entitled "Adaptation of linearized equation for vortex motion to 
twin -vortex model," the unperturbed fields v° and H° at one vortex will be those due 
to the other vortex. For vortex A in figure 11, equation (67) becomes 



where the brackets have been dropped as in equations (49) and (56). And for vortex B, 


P o (u B - v B ) 


x f B + p 0 h b x T b = o 


(68b) 


If the cylindrical coordinate system R, \p, z, with z perpendicular to the ground (see 
appendix A), is applied to both hydromagnetic vortices in figure 11, the radial components 
of equations (68) become at the ground 


Po( u ^-^) I t + ^ H X =0 

Po( U 5 - + '•’‘K'z = 0 


(69a) 

(69b) 


Formulas for and at the ground plane near a tilted hydromagnetic vortex 

are derived in appendix A. The fluid velocity fields v^; and v B at the ground plane 
are found by applying equation (A24) of appendix A 


2ttR c 

where R,, is the center -to -center distance in the ground plane and 

-p 

>'B = -#( 1 - sin X COS, M 

D = \Jl - sin^xcos^i^ ^1 + tan^x sin^i p^ 

Also, 

B r A 

^ 27TRC 

r A 

y A = — (1 + sin y cos i // A ) 


(70a) 


(70b) 

(70c) 


(70d) 

(70e) 
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A R 

The magnetic fields and are obtained from equation (A2 3): 


rx A _ X B 
H ^"2^ 

where 

I B 

= D^ 1 ' Sin X C ° S * A ) 

and 

H B = *A 
^ 27TRC 

where 

I. 

i A = -^(1 + sin x cos i^ A ) (7 Id) 

Expressions (70) and (71) will be used as mean values for the velocity and magnetic fields 
over the vortex cross sections. 

The vortex end points A and B are assumed to revolve about the point P (fig. 11) 
defined by the equations 


(71a) 


(71b) 


(71c) 


y A r A = ^B r B 


(72) 


r A+ r B = Rc 


with displacement velocities 



= r 




A r A 


(73) 


(74a) 


and 


< 14b) 

where the revolution rate ^(1) = ^gCO is now a function of the axial electric currents 
in the vortex cores. Equation (72) may also be written, by use of equation (73), 

(75) 

r A «c 
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Solution for revolution rate of inclined hydromagnetic v ortex pair at the ground . - 
Substituting equations (70), (71), and (73) in equations (69) and noting that T z = r cos y 
and I z = I cos x for both touchdown points (see appendix A) gives at A 


and at B 


'o r A^-'A*A»] = 
P o r B[d^- r B lp B (I *] ■ 


^A 1 ] 

27tR« 


^o I B i A 

2wRc 


(76a) 


(76b) 


If equation (76a) is multiplied by the 
by the factor (l - sin y cos 
tions (70b), (70e), (71b), and (7 Id) 


factor (1 + sin y cos and equation (76b) 

equations (76) become upon substitution of equa- 


r r B „ /tH _ ,i ° i A i B 

^A^ttRc r A^A^J 27711c 

and 

P o y B|p^ - r B^B (I) J = 


(77a) 


(77b) 


Equating the left-hand sides of equations (77a) and (77b) and collecting terms gives 


y A r A^A a) = ^B^B® 


(78) 


Relation (72) then shows that 


!// b (I) = tf/ A ( I) (79) 

which verifies that point P is the center point of revolution of points A and B in the 
ground plane, as assumed for equations (74). 

If equation (77a) is multiplied by 27 tRc and equation (75) is substituted, then 


p o y A y B 


2ttR^ a (I) 

y A +r B 


“ ^ i o 1 A 1 B 


(80) 
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But 


is the revolution rate ^ of touchdown points A and B in the absence of 


2»R 1 

axial electric currents, and equation (80) may be written 


p o y A y B 


1 - 


’/'a a) 




= Wb 


( 81 ) 


where denotes the revolution rate of points A and B with axial electric currents, 

and the revolution rate without electric currents. Substitution of definitions (70b), 
(70e), (71b), and (71d) also gives 


p o r A r B 


PAD 

i - -r 


Pi 


~ /i o I A I B 


(82) 


Now take the axial electric current for each funnel to be proportional to its individ- 
ual circulation 


X A “ « r A 


X B ~ * r B 


(83) 


Substitution of these relations in equation (82) results in the following equation for k, 
in A-sec/m 3 , 


K =1 


PAD 


Pt 


(84) 


Order of mag nit ude of axial e lectric current which would account for retarded 
revolution rate of tw in tornado . - In appendix C it is estimated that the revolution rate of 
the twin funnels was retarded by at least one-half with respect to the fluid flow and that 
the actual circulation of the pair + Tg was probably greater than 10^ m 3 /sec. 
Substituting the values 


*A + I’b = m 3 /sec 


(85a) 


p Q = 1.18 kg/m 3 


(85b) 
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(jLq = 4n x 10~^ H/m 


(85c) 




*1 


1 

2 


(85d) 


in equation (84) gives for air near sea level 
k = 6.86 x 10^ amp-sec/m^ 


( 86 ) 


and by equations (83) the total axial electric current which would be required to slow down 
the free revolution rate of the twin tornado to one -half that of the fluid flow is 


! a + ! b = 6.86 x 106 A 


(87) 


It was estimated in the section entitled "Effect of Core Buoyancy on Revolution Rate 
of a Twin Tornado” that buoyant attraction could have reduced the revolution rate of the 
twin tornado by about 5 percent. To accomplish the same 5-percent reduction by mag- 
netic attraction would require a total axial electric current of 2.16 x 10 6 amperes in the 


com- 


funnels. This would make the total magnetic energy of a tornado §§§ ^ 

parable in magnitude with its kinetic energy ^ — dV. Such electric currents are 


much larger than values previously proposed as possibly existent in tornadoes (refs. 20 
to 24). 


CONCLUDING REMARKS 

The new equation for the motion of vorticity derived in this paper is found to have 
application to the motion of tornadoes and tornado cyclones. Quantitatively, the most 
successful of these applications of the equation in its linearized form appears to be the 
explanation given for the tendency of severe local storms (tornado cyclones) to move to 
the right of the mean tropospheric wind by virtue of their interaction with the jet stream. 
Also successful was finding that the effect of the boundary layer at the ground on the twin 
tornado of April 11, 1965, near Elkhart, Indiana, was to make the two funnels separate, 
as they were observed to do. The increasing separation of the twin tornado also appears 
to constitute evidence against the possibility of convergent flow in the funnels of the twin 
tornado except within the boundary layer at the ground. 


43 



Less successful, quantitatively, were attempts at explaining why the funnels of the 
twin tornado apparently revolved about each other at a much slower rate than the fluid 
flow. The retardation in revolution rate was estimated to be about 50 percent with 
respect to the fluid. Two mechanisms were found which could retard the revolution rate: 
(1) buoyant attraction of each partially rarefied tornado core in the radial pressure gra- 
dient of the other core, and (2) magnetostatic attraction between axial dc electric currents 
of like sense in the cores. The retardation due to buoyant attraction was calculated to be 
about 5 percent for the twin tornado. In general buoyant attraction could retard the revo- 
lution rate of hollow -core or helium-core vortices of like sense by as much as 22 percent. 
Retardation by magnetostatic attraction requires electric currents in the tornado cores 
which are much larger than existing evidence would indicate. 

A potentially important buoyancy effect was not included in this report, that is the 
upward lift on buoyant atmospheric vortices which are tilted with respect to the vertical. 
The most significant new result appears to be the equation for vortex motion, itself. 

Since it is free of constraints on the fluid medium or on the flow or force fields, it should 
find applications in meteorology concerning the motions of various cyclones. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 24, 1970. 
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APPENDIX A 


DERIVATION OF MAGNETIC AND VELOCITY FIELDS AT THE GROUND 
PLANE NEAR A TILTED HYDROMAGNETIC VORTEX 

Determining the velocity field due to a tilted vortex above the ground is identical 
with the problem of determining the magnetic field due to an electric current in a tilted 
wire above the ground when it is specified that the magnetic field does not penetrate the 
ground. This is assured by making the ground plane a perfect electrical conductor, as 
shown in figure 12(a). Of these two analogs, the electrical lends itself simply to solution 


z 



(a) Original configuration with tilted 
current-carrying wire terminating 
in ground plane. 



(b) Equivalent image configuration 
for determining the magnetic 
field for z ^ 0. 


Figure 12.- Image method for magnetic fields due to steady electric current I in tilted 

wire above perfectly conducting ground plane. 


by the method of images because 


- s P e - 

div J = , whereas div w = 0. 

8t 


This distinction 


means that a vortex line cannot terminate at an interior point of the flow field, while an 
electric -current filament can terminate (by allowing accumulation of charge at the ter- 
mination point) and this is what is needed in the method of images. Therefore, the 
electric -current — magnetic -field problem will be solved initially by the method of images, 
and the circulation — velocity -field solutions will be determined afterward by analogy. 


In the image method the perfectly conducting ground plane of figure 12(a) is 
replaced by the image of the wire, as depicted in figure 12(b). The magnetic -field 
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APPENDIX A - Continued 


solution for this sharply bent wire is then obtained by superposition of the magnetic -field 
solutions for the two straight semi -infinite wire sections. The first concern, therefore, 
is to determine the magnetic field about a straight semi-infinite wire bearing a uniform 
current I, as illustrated in figure 13. 





Figure 13.- Semi-infinite wire in spherical coordinate system, 

— > 

showing uniform direct current I along the £-axis and 
charge accumulation Q e at the origin. 

A suitable representation for the current density in the wire is given by 


J ? = I6a)6(7])H(C) (Al) 

where 5 is the Dirac delta function and H(?) is the Heaviside unit step function. The 
continuity equation for current and charge densities 


9J C 

9C 


+ 



(A2) 


requires a buildup of negative charge density p g at the termination point of the wire, as 
represented by the expression 


p e = -It6(£)6(?7)5(C) (A3) 

(This negative charge will be canceled out upon superposition by an equivalent positive 
charge at the end of the image wire.) Maxwell's equations in a medium with negligible 
electric and magnetic polarization are 
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APPENDIX A - Continued 


curl H = e 0 E + J 

(A4a) 

div e 0 E = p e 

(A4b) 

curl E = 0 

(A4c) 

div H = 0 

(A4d) 


where H has been equated to zero because the free current density J = and the 
displacement current density e Q E are constants. 

It is apparent from the geometry in figure 13 and from equations (A4b) and (A4c) 
that the electric field E is purely radial, is spherically symmetric, and is given by 

E = REg 

where 


E 


R ~ “ 


It 

47re 0 R 2 


Hence, 



I 

4ire 0 R 2 


(A5) 


(A6) 


The integral form of equation (A4a) on a spherical surface S 


of radius R is 


(£ H • di = jj (e 0 E + j) • R dS (A7) 

If the closed contour line l is taken to be a circle perpendicular to the wire and if the 
area of integration S is that portion of the spherical surface enclosed by l, as depicted 
in figure 14, equation (A7) becomes 


27tRH^, sin 0 = 1 + 


f 2,e Q R 2 E 
J 0 ° 


Ej^ sin 0 d© 


(A8) 


Integrating the last term of equation (A8) and collecting terms gives for the axisymmetric 
azimuthal magnetic field 


TT _ _I (1 + cos 0) 

<P 4ttR sin 0 


(A9) 
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APPENDIX A - Continued 





Figure l4.~ Surface of integration S for equation (A7) 
and its bounding circuit line 2 on the spherical 
surface R = Constant in the coordinate system R,0,cp. 

Now take the straight semi-infinite wire to lie in the xz -plane of a Cartesian coor- 
dinate system, inclined at an angle y to the z-axis, as pictured in figure 15. The 
radius vector R will now be confined to the xy-plane. The amplitude of the magnetic 
field (eq. (A9)) at all points in the xy-plane (or Ri/'-plane) will now be determined. 

All that is necessary is to evaluate cos 6 and sin 9 of equation (A9) in terms of the 
angles y and if/. 

The formula for cos 9 from analytic geometry is 


cos 6 = AjA r + PjPr + 


(A10) 


where the direction cosines of I and of line 


R are given by 


Aj = sin y 

Pl = ° \ 

Vl = cos y J 

and 

A r = cos if/ ^ 

Vjl = sin if/ ^ 

0 


(All) 


(A12) 
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z 



Figure 15. - Configuration of figure 13 tilted and embedded in a Cartesian 
coordinate system x,y,z. The current-carrying wire is taken to lie 
in the xz-plane and the radius vector R in the xy- plane. 


Thus, 


and 


cos 9 = sin x cos 


sin 9 - y 1 - sin^x cos 

and the magnitude of the magnetic field at points in the xy -plane is given by 


(A13) 


(A14) 


tt _ I (1 + sin y cos 4/) 

<P 4ttR r — j IT- 

y 1 - sin^x cos* 4 / 


(z = 0) 


(A15) 


This magnetic field is oriented in the direction of r where r is normal to the 
plane formed by the vectors T and R, as shown in figure 16. 
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APPENDIX A - Continued 


2 



Figure l6.- Direction t of magnetic field at z = 0 

due to a semi-infinite tilted current-carrying wire. 

The horizontal component, denoted H^2 , is also shown. 

If XfUjV are the direction cosines of f, the condition that t be perpendicular to 
both T and R is expressed by 


XjX + VjV = 0 

(A16a) 

ArA + = o 

(A16b) 

A 2 + jU 2 + f 2 = 1 

(A16c) 


where the fact that /jj = v R = 0 has been used. Substituting v from equation (A16a) 
and n from equation (A16b) into equation (A16c) and solving for A 2 gives 


X 2 = 1 + 


X I . X R 

T 


-1 


‘l % 

or from equations (All) and (A12), 


(A17) 


(A18) 
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APPENDIX A — Continued 


In figures 12 and 16, the z -component of at points on the xy -plane will be 

nullified by the z -component of the image wire. But the horizontal component, which is 
proportional to sin a, will be doubled by the image field on the xy -plane. Therefore, 
the total magnetic field in the ground plane xy due to current in a straight tilted wire 
above the ground plane is horizontal, purely azimuthal, and given by 


= 211^ sin a (z = 0) 

In terms of the direction cosines of r, 


sin 2 o! = 1 - y 2 = X 2 + p 2 


and from equation (A16b) 



(A19) 


(A20) 


(A21) 


Substituting equations (A12) and (A18) gives after some reduction 


sin a 


-(■ 


+ sin^i// tan^x) 


- 1/2 


(A22) 


Substitution of equations (A15) and (A22) in relation (A19) gives for the magnetic 
field in the ground plane due to current I in a straight tilted wire above the 

ground plane 


H ^/r 

where the coordinates R 
in figure 17. 


1 + sin x cos i// ( z _ ( 

sin 2 x cos 2 ^ yl + tan^x sin 2 i// 

and if/ in the ground plane and the tilt angle 


X 


(A23) 


are shown 
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APPENDIX B 


DESCRIPTION AND TABULATION OF EMPIRICAL DATA 
ON TWIN TORNADO OF APRIL 11, 1965, 

NEAR ELKHART, INDIANA 

An unusual tornado, identified J-2 in reference 4, occurred on April 11, 1965, near 
Elkhart, Indiana. This tornado developed a second funnel, and the two funnels revolved 
about each other for about 1 minute, and then the original funnel disappeared. Six photo- 
graphs of the tornado were taken during the period of duality by Paul Huffman of the 
Elkhart Truth. One of these photographs is reproduced in figure 5. On the basis of these 
photographs and aerial examination of the damage paths, Fujita determined that the circu- 
lation of both funnels during their brief period of individuality was cyclonic, and he pro- 
duced figure 4 (fig. 49 in ref. 4) which shows to scale the motion of the two vortices. The 
time intervals between the photographs were determined from the translational speed of 
the storm, about 22.4 m/sec (50 mph). 

From a fourfold enlargement of figure 4(b), the following parameters have been 
obtained for each of the locations 3 to 6: (1) the center -to-center distance Rc of the 

twin vortices; (2) the tilt angle x °f the pair with respect to the vertical; and (3) the 
azimuth angle ip^ of the pair in the ground plane. These quantities are listed in 
table I. The aximuth angle xp^ is the angle at funnel A measured counterclockwise 
from the broken line labeled "tilt 29°" in figure 4(b) to the line AB connecting the cen- 
ters of the vortex pair. The diameters of the spots in figure 4(c) were taken to be indic- 
ative of their relative strengths and were used to calculate the tabulated values for the 
ratio r B /r A . 

For each of the three intervals between these four data points interpolation gives 
the mean separation (Rc)> the mean tilt angle (x)> the mean azimuth angle 

V r 

quantities are listed in table H. Also included for use in appendix C are mean values 
<D) where D = yl - sin^x cos^i/^^l + tan^x sin^^. 


A ^. These 


the mean revolution rate the mean circulation ratio ^ 
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APPENDIX C 


CIRCULATION OF THE TWIN TORNADO AS DEDUCED 
FROM ITS REVOLUTION RATE - 

EVIDENCE THAT THE FUNNELS MOVED MORE SLOWLY THAN THE FLUID 


Circulation Estimate From Revolution Rate 

In this appendix the total circulation + Tg of the twin tornado will first be 
estimated from its revolution rate with the assumption that the funnels move with the 
fluid. According to equation (36) this means that the buoyancy and source terms are 
neglected as well as a possible Lorentz force F. The effects of these terms have been 
determined, subject to certain simplifications, in the sections entitled, "Effect of Core 
Buoyancy on Revolution Rate of a Twin Tornado,” "Effect of the Boundary Layer at the 
Ground on Tornado Motion," and "Influence of Axial Electric Current on Revolution Rate 
of Twin Vortices." After the first estimate is obtained, the effects of the neglected terms 
on this estimate will be considered. 

The velocity fields are determined in accordance with the twin -vortex model pre- 
sented in the section entitled, "Tornado model" and illustrated in figure 8. Effects of 
the tilt angle x and of the unequal circulations and Tg of the two vortices are 

included. But the convergences and Qg of this model are equated to zero in order 
to conform to the observation that the twin funnels did not approach each other, as 
explained in the section entitled, "Correlation of the twin-vortex model with data from the 
twin tornado." 

Velocity field at ground plane due to single inclined vortex . - Consider a single vor- 
tex of circulation r which terminated on the ground plane and is inclined to the vertical 
by the angle x as in figure 18. The velocity field at the ground plane (Ri//-plane) for 
such a vortex in a perfect fluid is found in appendix A (eq. (A24)) to be purely azimuthal 
and given by the formula 


v^(R,^) = 


2ttR 


1 + sin y cos ^ 


\fi - sin^x cos^i/z/l + tan^x sin^ 


(z = 0) 


Twin -vortex model . - Now, consider the configuration shown in figure 8 as viewed 
from above in figure 19. Here the two tilted, semi -infinite vortices terminate at spots A 
and B in the ground plane. The broken lines in this figure represent the projections of 
the vortex axes in the ground plane, as in figure 4(b). Take the two inclined vortices to 
be parallel to each other and to have cyclonic circulations and Tg in a perfect 
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/ 



Figure 19.- View of ground plane from above (as in fig. Mb) ) 
showing the touchdown points of two vortices A and B, their 
velocity fields, and orientation with respect to projections 
of the vortex axes . 

fluid which is otherwise at rest. The velocity field v® at point B is perpendicular to 
the line AB, as shown in figure 19, and is given by 


B 7 A 

VTT = 


P 2irR n 


where 


7 A = "d ( X + sin ^ COS ^A.) 


and 


D = ^1 - sin 2 x cos 2 i^ a ^1 + tan 2 y sin 2 ^ 


Similarly, the velocity field v^ is given by 


T*. rB 


P 2nR n 


where 


7 b = iri 1 - sin * cos *a) 


(Cl) 


(C2) 


(C3) 


(C4) 


(C5) 
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and the substitution = ^ A ± it has been made. The quantities y A and y^ will be 
referred to as the "effective circulations" of vortices A and B. 

Define point P on the line AB, as shown in figure 19 so that 


^A r A =r B r B < C6a > 

r A + r B = R c (C6b) 


Then the revolution rates of points A and B about P are given by 


^A 



^B 

27rR c r A 


^B 



y A 

27rR c r B 


But by relation (C6a) 


^B 


= ^A 


(C7) 


(C8) 


(C9) 


and point P is seen to be the center point of revolution of points A and B in the ground 
plane, (since y A and y B are different functions of s// A , the location of point P on 
line AB will also vary with the azimuth angle of the pair.j Combination of equa- 
tions (C6a) and (C6b) gives 


^B _ ^A + ^B 
r A R c 


(CIO) 


Hence, the theoretical revolution rate in the ground plane of a tilted vortex pair of like 
sense in a perfect fluid is given by 


r A + 

- J ^~r (C11 > 


The effective circulation y A + y B of the vortex pair is related to the circulation 
T a + Tg by adding equations (C2) and (C5) 
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y K + ^B = 


r A +r B 

D 



sin x cos i// A 


(C12) 


where D is given by equation (C3). Combining equations (Cll) and (C12) gives 


2ttR^D^ a = r A (l + b) + r A (l - b)sin x cos \p A 


(C13) 


where b = Solving for r A : 


r A = 


27rR^D^ A 


(1 + b) + (1 - b)sin x cos ^ A 


(C14) 


and multiplying by (1 + b) gives 

27tR 2 D^ a (1 + b) 


p + P - 

" B (1 + b) + (1 - b)sin x cos i// A 


(C15) 


For intervals 3-4, 4-5, and 5-6 mean values ^r A + can be determined by substi- 
tuting values from table II for (Rc)’ (D), <^ a ), (b), (x)> an d <^A> in 

formula (C15). This gives the following estimate for ^r A + Tg^ in mVsec: 


5.51 x 10 3 (Interval 3-4) 

< r A + r B) = S 4,34 X 103 ( Interval 4 ' 5 ) 

4.91 x 10 3 (Interval 5-6) 

v. 

Averaging these values gives 

( r A + r B> = 4,92 X 103 ( interval 3 " 6 ) 


(C16) 


(C17) 


for the total circulation in m 3 /sec of the tornado pair, as computed from the revolution 
rate of their touchdown spots on the ground with the assumption that the funnels move 
with the fluid. 
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Evidence That the Revolution Rate Was Retarded 


With Respect to the Fluid 


Examining the wind damage shown in figure 20 shows that the value (C17) for 
(r A + r B ) is too small by at least a factor of 2. This photograph views the Midway 
Trailer Court from the north. This trailer court was hit by the twin tornado. The two 
funnels moved from right to left across this photograph and were tilted toward the bottom 
edge of the photograph. Figure 5, which depicts the funnels at location 4, was taken while 
funnel B was over the trailer court. 



Figure 20.- Aerial photograph of severe damage to the Midway Trailer Court near 
Elkhart, Indiana. Devastation seen near center of picture was caused by the 
left funnel of the tornado shown in figure 5> while the right funnel moved 
over the plowed field from right to left, cutting across upper comer of the 
court. Damage occurred at 5:32 p.m. CST, April 11, 1965; picture was taken 
April 13. Courtesy T. Fujita (ref. 4). 


Note the damage near the lower edge of figure 20 left of center. One trailer has a 
side blown out, and next to it is one which was upturned (wheels up). The maximum pos- 
sible wind speed at this point due to the calculated value (C17) for ^r A + can be 

determined by combining the two funnels into a single vortex on the nearest tornado path 
and tilting this vortex directly toward the upturned trailer. The center of the nearest 
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path is a distance of 89.9 m (295 ft) from the upturned trailer. (This distance estimate 
assumes that the trailers were typically 15.2 m (50 ft) long.) The nearest path is that of 
funnel B, which cut diagonally from upper right to center left across the trailer court and 
caused the devastation near the middle of the photograph. 

The velocity field at the ground due to a single tilted vortex (fig. 18) is given by 
equation (A24). Substituting the values 

r = (r A + r B ) = 4.92 x 10 3 m2/sec 
R = 89.9 m (295 ft) 

X = 29° 

in equation (A24) gives 

v^= 14.8 m/sec (33.1 mph) (C19) 

as the maximum possible wind speed at the upturned trailer due to the circulation value 
< r A + r B> ' 4 92 x 10 3 m2/sec calculated from the revolution rate of the pair, assumed 
to move with the fluid. This wind speed is clearly too low, and it is not augmented by the 
speed of the storm, about 22.4 m/sec (50 mph), because the site of the upturned trailer 
was on the northern side of the tornado path, where the storm velocity opposes and tends 
to reduce the circulation velocity. In reference 3, Fujita reports an anemometer 
recording which shows winds of 67 m/sec (150 mph) at the edge of a damage path of a 
tornado generated an hour later by the same storm. The conclusion therefore is that the 
actual circulation of the twin tornado was much greater than the value (4.92 x 10^ m2/sec) 
calculated from its revolution rate by assuming that the vortices moved with the fluid. 
Apparently, the revolution rate was retarded with respect to the fluid by at least a factor 
of 1/2. 

The possibility of local wind shear causing this retardation has been considered. 

A continuous source of wind shear is the tornado cyclone in which tornado J-2 was 
embedded. This wind shear can be estimated from the paths of the family of tornadoes 
generated by the tornado cyclone, as presented in the maps of reference 4. The wind 
shear estimated in this way is found to be too low to affect the revolution rate of the twin 
tornado J-2 significantly. 
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TABLE I.- PARAMETERS FOR TWIN TORNADO J-2 
OBTAINED FROM FIGURE 4 


Location 

Center -t o -c enter 
separation in 
ground plane, 

*c 

Tilt angle with 
respect to vertical, 
X , de S 

Azimuth 
angle, 
^ A , deg 

r B / r A 

ft 

m 

3 

390 

118.9 

31 

96.9 

0.6 

4 

431.8 

131.6 

29 

134.1 

1.00 

5 

478.2 

145.7 

29 

151.4 

1.35 

6 

482.6 

147.1 

29 

237.0 

2.75 


TABLE II.- MEAN VALUES FOR PARAMETERS OF TABLE I AND 
MEAN OBSERVED REVOLUTION RATES FOR TORNADO J-2 


Interval 

between 

locations 

Mean 

separation, 

<«c> 

Mean tilt 
angle, 
<X>, deg 

Mean azimuth 
angle, 
<^ A ) , deg 

Mean revolution 
rate (observed), 
<^A> ’ r a d/sec 

< r B/ r A> 

<D> 

(eq. (70c)) 

ft 

m 



3-4 

(13 sec) 

410.9 

125.2 

30 

115.5 

0.05 

0.8 

1.1 

4-5 
(8 sec) 

455.0 

138.7 

29 

142.8 

.038 

1.17 

.973 

5-6 

(31 sec) 

480.4 

146.4 

29 

194.2 

.048 

2.05 

.891 
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